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ABSTRACT 

We present a new distance determination to the Galactic globular cluster 47 Tucanae by fitting the 
spectral energy distributions of its white dwarfs to pure hydrogen atmosphere white dwarf models. 
Our photometric dataset is obtained from a 121 orbit Hubble Space Telescope program using the Wide 
Field Camera 3 UVIS/IR channels, capturing F390W, F606W, F110W, and F160W images. These 
images cover more than 60 arcmin 2 and extend over a radial range of 5 — 13.7 arcmin (6.5 — 17.9 pc) 
within the globular cluster. Using a likelihood analysis, we obtain a best fitting unreddened distance 
modulus of (m - M) Q = 13.36 ±0.02 ±0.06 corresponding to a distance of 4.70 ±0.04 ±0.13 kpc, where 
the first error is random and the second is systematic. We also search the white dwarf photometry 
for infrared excess in the F160W filter, indicative of debris disks or low mass companions, and find 
no convincing cases within our sample. 

Subject headings: Galaxy: globular clusters: individual(47 Tucanae) — white dwarfs — stars: dis- 
tances 



1. INTRODUCTION 

Many astrophysical analyses involving globular clus- 
ters rely on well determined metallicity, reddening, and 
distance measurements. In particular, the estimated ages 
of Galactic globular clusters can depend heavily on these 
assumed quantities, so their careful determination is nec- 
essarily important. Ages of Galactic globular clusters can 
place strong constraints on the epoch of formation of our 
Galaxy, identify when cluster formation occurred with 
respect to the era of reionization in the Universe and the 
cosmic star formation peak, constrain the formation his- 
tory of the Gala xy by dating the halo and bulge, as well 
as constrain VLm (jGratton et al.ll2b03[ ) when coupled with 
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measurements of H by WMAP (iSpergel et ai1l2003D and 
the HST Key Project (jFreedman et al.ll200lD . Most of 
the Galactic globular clusters have ages between 9- 14 Gyr 
(e.g. IMarin-Franch et al"1l2009t IDotter et al.ll2010D . with 
most absolute ages commonly determined by fitting the- 
oretical isochrones to the main sequence turnoff. This, 
however, depends strongly on the assumed distance mod- 
ulus to the cluster, which is t he most significant co ntrib- 
utor to the age uncertainty (jBolte fc Hoganl [19951 ) . Al- 
though less sensitive to distance, even the more recently 
employed age dating technique via the white dwarf cool- 
ing sequence (IHansen et al.l F2004. 2007) requires a precise 
distance modulus. 

The Galactic globular cluster 47 Tucanae (47 Tuc, 
NGC 104) is considered the metal-rich proto-type for 
globular cluster wor k ([Fe/H] = -0.74 ± 0.003 ± 0.026, 
iCarretta et al.ll2009l ). As the closest and least reddened 
metal-rich globular cluster, 47 Tuc serves as an impor- 
tant science, as well as a calibration target, for the 
Hubble Space Telescope (HST). It is used as a metal- 
rich anchor when comparing observations of resolved 
stars in nearby galaxies to stellar evolutionary models 
(iCassisi fc Salarisll997l:lZoccali et al.lll999tlFerraro et al.l 
119991 IDi Cecco et al.l I2010D and in examining the star 
formation history o f galaxies in the Local Group (e.g. 
IMonelli et al.l 120101) . As another example, the Wide 
Field Camera 3 (WFC3) Galactic Bulge Treasury pro- 
gram (GO-11664) has ta rgeted 47 Tuc in t heir resolved 
stellar population study ()Brown et al.ll2009l ). With these 
observations, ranging from the UV to near-IR, 47 Tuc 
will serve as an empirical population template, helping 
to calibrate the reddening-free indices and correcting the 
transformation of theoretical isochrone libraries into the 
WFC3 photometric system. 

With relatively well measured values of metallicity 
and reddening (E( B - V) = 0.04 ± 0.02, iSalaris et all 
l2007t lHarrisI 12010). the distance becomes the crucial 
input parameter to determine the age of 47 Tuc. A 
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TABLE 1 

Comparative Absolute Distance Estimates to 47 Tuc 



(m - M) a 



Method 



Reference 



13.36 ± 0.02 ± 0.06 b WD Spectral Energy Distributions 

13.23 ± 0.08 Eclipsing Binary 

13.25 Main Sequence Fitting 

13.32 ± 0.09 Tip of the Red Giant Branch 

13.47 ±0.11 RRLyrae 

13.28 ± 0.07 Eclipsing Binary 

13.18 ± 0.03 ± 0.04 b Horizontal Branch Fitting 

13.02 ± 0.19 Cluster Kinematics 

13.40 ± 0.08 c ' d Main Sequence Fitting 

13.25 + 0.06 - 0.07 Main Sequence Fitting 

13.21 ± 0.04 ± 0.07 b Main Sequence Fitting 

13.15 ± 0.14 WD Cooling Sequence Fitting 

13.45 ± 0.09 c Main Sequence Fitting 

13.32 ± 0.20 Tip of the Red Giant Branch 

13.29 ± 0.20 c Horizontal Branch Fitting 
13.57 ± 0.15 Main Sequence Fitting 
13.32 ± 0.03 ± 0.04 b Horizontal Branch Fitting 
13.23 ± 0.17 Baade-Wesselink of RR Lyrae 
13.38 ± 0.05 Horizontal Branch Fitting 
13.52 ± 0.08 c Main Sequence Fitting 
13.28 Horizontal Branch Fitting 
13.18 Main Sequence Fitting 



This Paper 

Thompson et al. (2010) 
Bergbusch fe Stetson (2009) 
Bono et al. (2008) 
Bono et al. (2008) 
Kaluznv et al*T 2007) 
Salaris et al. (2007') 
McLaughlin et al. (20061 
Gratton et al. (2003) 
Percival et al. (2002) 
Grundahl et al. (2002) 
Zoccali et al. (2001) 
Carretta et al. (2000) 
Ferraro et ld"72000) 
Ferraro et al. (1999) 
Reid (1998)~ 
Kaluznv et al. (1998) 
Storm et al"7l994) 
Salaris fc W eiss (1998) 
Gratton et al. (1997) 
Hesser et aT*Tl987r 
Hesser et al. (1987) 



a Dereddened using E(B - V) = 0.04 and A v / E(B - V) = 3.1. 

b The first uncertainty is random and the second is systematic. 

c Derived from the quoted distance modulus uncorrected for binary bias. A binary 

correction removes 0.02 mag from the distance modulus. 

d Using their mean distance derived from two colors. 

e From their global metallicity derived distance modulus. 



considerable effort has gone into determining the dis- 
tance to this cluster. We list recent work in Ta- 
ble [TJ There has been quite a variety of method- 
ologies used to determine the distance of this cms 



ter, including main s equen c e fitting dHesser et ail 
Gratton et aU Il997t iReidl 119981: ICarretta et"a!T 



1987 



2000; 



Grundahl et al.ll2002l; rPercival et al.ll2002l: IGratton et all 
20031 : [Bergbusch & Stetson 2009) , the tip of t he red giant 
branch (iFerraro et alll2000t [ Bono et al. 2008), horizontal 
branch fitting (IHesser et aLl ll987t iKaluzny et all 119981: 



Salaris fe W eiss! 119981: IFerraro et all 119991: ISalaris et all 



20071), RR Lvrae (IStorm et allll994HBono et all | 2008[) 



cluster stellar k inematics ([McLaughlin et al. 12006), 



eclips ing binaries ([Kaluznv et all 120071 : iThompson et all 
I2010D. and white d warf (WD) cooling sequence fitting 
(|Zoccali et al.l [200lT) . This set of measurements spans 
unreddened distance moduli between 13.02 — 13.57 for 
47 Tuc which has important implications for the derived 
age. Examining only the main sequence turnoff lumi- 
nosity for a cluster at the extremes in distance moduli, 
this yields a difference in age of ~ 6 GviPl, altering any 
inferred history of the formation of the Milky Way that 
one might derive from 47 Tuc. 

Our goal is to determ ine the distan c e to 4 7 Tuc fol- 
lowing the example of iZoccali et all (120011) by using 
the cluster's WD population. IZoccali et alT(|2001l ) con- 
structed a WD sequence using local calibrators with mea- 
sured trigonometric parallaxes and magnitudes to com- 

11 To first order, we find that when examining the turn off 
magnitude, the age uncertainty is approximatley equal to (10 
Gyr /mag)* (uncertainty in the distance modulus). 



pare with the WD cooling sequence of 47 Tuc, yielding an 
unreddened distance modulus of (m - M) D = 13.15 ±0.14 
(see Tabled]). Here, rather than use local calibrators, we 
determine the distance by comparing our measured pho- 
tometry to the hydrog en atmosphere (DA) WD models of 
iTremblav et all (12011) and the helium atmosphere (DB) 
models of lBergeron et al.l (|201lD . Previous spectroscopic 
studies have found solely DA WDs in globular clusters 
within the temperature r ange of our samp le, including 
the globular clus ters M4 (iDavis et al l 120091) . NGC 6752 
and NGC 6397 dMoehler et al.l 120041) as w ell as in the 
open cluster NGC 2099 (jKalirai et alJl2005T ). While this 
has not been tested in 47 Tuc to date, we make the as- 
sumption that our sample consists of DA WDs, but we 
also compare our results to that obtained considering DB 
models and a mixture of DA and DB WDs. 

WDs within globular clusters are faint objects in 
crowded fields, making their identification quite chal- 
lenging. However, with carefully planned observations 
with the HST, combined with thorough data reduction 
and analysis, it becomes possible to obtain high quality 
multi-band photometry of WDs in large samples within 
one cluster system (see t he recent observations of WDs in 
M4 and NG C 6397 fromlRicher et al.l l2002HHansen et all 
l2lMl2lMlRicher et alll2008D . 



2. DATASET 

We obtained HST images of the Galactic globular clus- 
ter 47 Tuc in cycle 17 (GO-11677, PI - H. Richer) with 
the Advanced Camera for Surveys (ACS) and the WFC3. 
The primary field of this program was imaged with the 
ACS/WFC in F606W and F814W filters for 121 orbits. 
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Additional fields were taken with the WFC3 using both 
the UVIS and IR channels. Our obs ervational des i gn an d 
reduction are discussed in detail in lKalirai et al.l (|2011l) . 
here we briefly review the WFC3 dataset used in the 
present analysis. With the WFC3, we obtained 121 or- 
bits covering > 60 arcmin 2 and extending over 250 de- 
grees in azimuthal range of 47 Tuc. This covers a ra- 
dial range of 5 — 13.7 arcmin (6.5 — 17.9 pc) within the 
globular cluster. For 61 orbits, we stared in one field lo- 
cation in the cluster (stare field) obtaining 59 orbits of 
IR images in F110W and F160W and 2 orbits of UVIS 
images in F390W and F606W. We then swept out an 
arc around the cluster's center (swath fields), using 5 or- 
bits in each of our 12 additional fields obtaining F390W, 
F606W, F110W and F160W image s (for a visualizatio n 
of our observations, see Figure 4 of IKalirai etll|[20rl . 
The observational logs of t he swath and stare fi elds can 
be found in Tables 2 & 3 of IKalirai etal! (|2011D . 

To summ a rize the data analysis discussed in 
IKalirai et al.1 (l20Tl . we retrieved the WFC3 observa- 
tions from the MAST archive and re-processed all data 
into _flt images using the latest calibration files. These 
images were then supersampled and drizzled together, 
paying careful attention to image registration, sky offsets, 
and pixel drop size. PSF-fitted photometry and astrome- 
try was performed on the final image stacks in each filter 
using DAOPHOT II and ALLSTAR. The IR swath data 
were too undersamplcd (with a native FWHM= 1.1 — 1.2 
pixels) for PSF photometry however, so photometry for 
each star was obtained from a small aperture radius of 
3.5 pixels. The catalogs from all filters, over both the 
UVIS and IR cameras, were merged together using high 
order transformations. 

3. SELECTING THE SAMPLE 

In the selection of our sample of WDs, our goal was 
not to have the largest sample possible with the dataset, 
but instead to have a set of WDs that cleanly lay on 
the white dwarf cooling sequence and are free of possible 
photometric contamination. Therefore, we selected our 
sample using multiple criteria. Initially, we isolated the 
WD cooling sequence in the (F390W-F606W, F606W) 
CMD shown in Figure [TJ selecting WD candidates down 
to F606W = 25.5. 

Our WD cooling sequence was carefully selected to 
avoid possible contaminates from the main sequence of 
the SMC in the background. From our inital list, there 
were 171 WD candidates that also had F110W and 
F160W measurements. From these 171 candidates, we 
selected objects that were consistent with being stellar 
and passed a visual inspection of the images. 

To classify an object as stellar, full width half maxima 
(FWHMs) were measured in all 4 images and compared 
to the mean stellar FWHM in these same images. To de- 
termine the mean stellar FWHM in each image, we mea- 
sured the FWHM of 20-30 stars that lay on the main 
sequence between F606W = 19 - 19.5 in the (F390W- 
F606W, F606W) CMD. We selected the same 20-30 stars 
for the 4 filters in each field. We removed candidate ob- 
jects from our list of 171 candidates that differed in their 
measured FWHM by more than 1 standard deviation 
from the mean stellar value in the F160W filter. We 
chose the F160W filter because the resolution of the im- 
ages are the lowest in our sample and the WDs are thus 




F390W-F606W (mag) F390W-F60** (meg) 

Fig. 1.— Left: The UVIS CMD of the objects in the stare and 
swath fields. The reddest sequence is the main sequence of 47 Tuc 
and the bluest sequence is the cluster's WD cooling sequence. The 
intermediate sequence is the main sequence and turn off region of 
the SMC in the background. Right: A closer examination of the 47 
Tuc WD cooling sequence. Our final selection of 59 WDs are shown 
as green (light grey) circles while the remaining 112 candidate WDs 
are shown as red ( dark grey ) circles. A color version of this figure 
is available in the electronic version. 



more likely to be contaminated by nearby neighbors. We 
are also interested in searching for IR excess in the WD 
photometry. This may be indicative of a debris disk or 
a low mass companion (see Section RTTj) . While most of 
the IR excess, if present, will be evident at wavelengths 
longer than F160W, this program could be able to see 
hints of the excess in the F160W band which may be 
followed up in future observations. As such, we are par- 
ticularly critical of the F160W measurements. The dis- 
tributions of stellar FWHMs and the final selected WDs 
are shown in Figure [2] 

Independent of our size rejection analysis, we also visu- 
ally examined the F390W, F606W, F110W, and F160W 
images for the 171 WD candidates to search for objects 
that appeared elliptical or whose photometry may have 
been contaminated by bright neighboring stars, diffrac- 
tion spikes, remaining cosmic rays, or those found in ex- 
tremely crowded regions with varying background. If the 
candidate object was contaminated in any of the four fil- 
ters, it was rejected from our list. 

We found in these two search methods that all of our 
visually rejected objects were also rejected using our 
FWHM technique. For a few objects, we were unable 
to measure a reliable FWHM in F160W, so we accepted 
these objects into our final list if their FWHMs were 
stellar in the three additional bands and they did not 
appear contaminated on the images. Our final selection 
consists of 59 WDs listed in Table [2] with their field lo- 
cation, R.A. and Deck (in J2000 epoch), their F390W, 
F606W, F110W, and F160W photometric measurements 
and uncertainties in VEGAMAGS, and the effective tem- 
peratures that we fit in this study. In Figure EH we show 
the F390W, F606W, F110W, and F160W images for our 
sample WDs. Two of these WDs were found near the 
edge of their F390W and F606W images. For these two 
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Fig. 2. — The measured FWHMs for main sequence stars between 
F606W = 19 — 19.5 magnitudes in all fields are shown as the open 
histogram for the filters (top left) F390W, (top right) F606W, (bot- 
tom left) F110W, and (bottom right) F160W. The solid histograms 
are the measured FWHMs for the 59 selected WDs. 




Fig. 3. — A four panel image of WD 1. Within the four panels, the 
upper left is F390W, upper right is F606W, lower left is F110W, 
and lower right is F160W. The WD is centered in the 2" x 2" 
images and its photometry is listed in Table [2] 

filters, there were only two exposures combined to pro- 
duce the final image and thus some cosmic rays were not 
removed due to the offset between the two exposures. 
This is evident in Fig. [3J for WDs 23 and 52 but does not 
affect our photometry in these bands obtained with PSF 
fitting. 

4. SPECTRAL ENERGY DISTRIBUTIONS 

4.1. Model Conversion to our Photometric System 

In order to generate the spectral energ y distributions of 
WDs i n 47 Tuc, we have chosen to use the lTremblav et afl 
(|2011f ) spectral model grids of pure hydrogen atmo- 
sphere WDs. These model grids define spectra for ob- 
jects with surface gravities (log(g)) ranging from 6 — 10 
in steps of 0.5 and effective temperatures (T e gp) from 



6000 - 120000 K. The spectral grid is defined in steps 
of 500 K between 6000 - 17000 K, steps of 5000 K be- 
tween 20000 - 90000 K, and steps of 10000 K between 
90000 - 120000 K. We have converted the model grids 
to our observed Space Telesc ope magnitudes (STMAG), 
similar to the procedure of IBergeron et ail ([1997D . to 
obtain a magnitude for each of our observed filters, 
•MstmAG- This enables us to fit the models to our obser- 
vations via Equation [TJ 



M s 



TMAG 



-2.5 log 



IF 



10 



fo° ^FmodE\S\ dX 



J™XF S x dX 



(1) 



where F mo d is the flux per unit wavelength for the spec- 
tral model integrated over the solid angle, 5a is the total 
system throughput for the filteiF l. and E\ is the stan- 
dard interstellar extinction curve P?itzpatricklll999|) with 
A v I E(B - V) = 3.1, using some input E(B - V). The 
STMAG zeropoints are defined so that a flat spectrum 

of 3.63xl0 -9 ergs cm~ 2 s _1 A (Frj) through any filter 
yields a magnitude of zero. 

The inclusion of the A terms in Equation [T] is required 
because of how the flux is measured on a CCD. The CCD 
does not know the individual energy of each photon that 
it receives, only whether or not it has received a photon. 
For any given energy density, the number of photons re- 
ceived will increase linearly with wavelength. For a rela- 
tively flat spectrum, fluxes measured on a CCD through 
some filter will therefore be more affected by the energy 
flux at longer wavelengths. This is because this region 
of the spectrum will necessarily have a higher number of 
photons per unit energy as the energy of a photon is in- 
versely proportional to its wavelength (E p hoton = ch/A). 

The factor of (R/d) 2 in Equation [T] scales the model 
flux, the flux per unit wavelength at the surface of the 
object, to the observed flux for a star of radius R, at a 
distance d. Equation Q] can then be rewritten as, 



Mstmag = —2.5 log 10 | 



J™XF S x dX j 
+5l0gl °(^ 



(2) 



Equation [21 which is illustrated in Figure HI shows that 
the model magnitude in any given filter is a function 
of only log(g), T e ff, reddening, and d/R. We apply this 
method to convert the model grid magnitudes to our four 
filters (F390W, F606W, F110W, and F160W). 

4.2. Fitting the Model SEDs to the WD Photometry 

Using WD cooling models from iFontaine et al.l ([2001), 
we can derive a relationship between surface gravity and 
effective temperature for a given mass (see Figure [5j. 
Additionally, for any given mass and log(g), we can cal- 
culate the radius from Equation [31 



R = 



IGM 



(3) 



12 The STMAG zeropoints and throughput values for each 
filter have been obtained from the Space Telescope Science 
Institute website l |http: //www.stsci.edu/hst/wfc3/phot_zp_lbn| 
http: / / www.stsci . edu / hst / wtc3 / ins_pertormance / throughputs) . 
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Fig. 5. — Left: Contours of constant surface gravity (log(g)) in 
mass-effective temperature (T fj) space and Right: Contours of 
cons tant radius, measured i n units of solar radii, in mass-T c g space 
from lTremblav et al.l pOTH V 



Fig. 4. — An illust ration of Eq uation [2] for a sample model spectrum converted to an F606W magnitude. Upper panel: An example 
model spectrum from Tremblay ct al. (2011) (F moc {) for a DA WD with an effective temperature of 10 000 K and surface gravity of 7.50. 
Overlaid are the four filter throughputs S\ (F390W, F606W, F110W, and F160W). Lower left: The interstellar extinction. Lower right: 
The combined spectum. 

Putting the mass dependent T c ff - log(g) relation to- 
gether with Equation [3] means that, of the four param- 
eters, T e ff, log(g), mass, and radius, the combination of 
T c ff and any of the other three allows one to calculate 
the remaining two parameters. Of these four, we have 
chosen to parametrize our grid over T c fj and mass, since 
previous work allows us to put a reasonable prior on the 
mass of our objects. This reduces the six parameters 
that control each WD spectrum (effective temperature, 
mass, distance, reddening, surface gravity, and the ra- 
dius of the WD) to four parameters that need to be fit 
in our procedure. We are therefore able to generate a 
model SED in these four filters given a mass, effective 
temperature, reddening, and a distance. The first three 
parameters control the shape of the SED, while the dis- 
tance controls the overall scaling. In order to determine 
temperatures for our WDs and to display both model 
and WD SEDs, we have first held mass and reddening 
as constants in our fitting procedure. We have assumed 
a mass of 0.53 M^. Thi s follow s fr om the the prediction 
of iRenzini fc Fusi Peccil (|1988f ) and lRenzini et all (119961 ) 
that the WDs forming today in globular clusters should 
have a mass of 0.53 ± 0.02 M Q . Further support of our 
mass selection com es from recen t direc t spectroscopic 
measurements from iKalirai et all ([2009) of 0.53 ± 0.01 
M for the mass of WDs in M4 and 0.53 ± 0.03 for 
the mass of WDs in NGC 6752 (|Moehler et alll2004l) . 
We also assume a reddening of 0.04 and then fit the dis- 



tance modulus and the effective temperature to each in- 
dividual object. We ignore the individual distance mod- 
uli here and include the fitted effective temperatures in 
Table H for our 59 WDs. Figure © shows the SEDs for 
our selected WD sample. 

In Scction[5j we determine the distance to 47 Tuc using 
Equation [5] to fit to our WD photometry, however we do 
so using a maximum likelihood function and marginalise 
over various parameters. 

5. DISTANCE TO 47 TUC 

To determine the distance to 47 Tuc, we will use white 
dwarf spectral models in order to predict an expected 
absolute magnitude for our sample WDs in our 4 filters, 
similarly to our technique in Section [4] However in this 
section, we fit the distance modulus to 47 Tuc by com- 
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5000 10000 15000 
Wavelength (A) 

Fig. 6. — The SED for WD 1. The filled squares are our pho- 
tometric measurements in F390W, F606W, F110W, and F160W 
with uncertainties. The wavelength range for each filter is also 
shown as the equivalent width of the filter divided by the maxi- 
mum throughput, as defined by the WFC3 instrument handbook. 
The red triangles, connected by the red line is the best fit model 
in STMAG magnitudes. A color version of this figure is available 
in the electronic version. 

paring the expected absolute magnitude to the measured 
apparent magnitudes of the WDs using a maximum like- 
lihood analysis that uses all 59 WDs to determine one 
final distance to the cluster. Unlike in Section 0] we do 
not assume a single value for reddening and WD mass 
in the determination below, rather we marginalise over 
each of these parameters with the inclusion of Gaussian 
priors. 



5.1. Maximum Likelihood Function 

Our procedure begins by calculating the likelihood 
function of our 59 WDs, termed Lj where j runs from 
1 — 59, over our four parameters of T er r, reddening, mass, 
and distance which is simply converted to a true distance 
modulus, (m - M) Q . This is shown in Equation^] 

Lj(data | (T cff , (m - M) ,mass,E(B - V))) = (4) 



n° xp 



(datai — mod, 



Y2l 



2a 



The photometric data for each WD (termed datai) are 
the measured magnitudes in our 4 filters, modi are 
the magnitudes calculated from the model described in 
Equation [2] from Section |4j and o~i are the uncertainties 
for each point. 

After a likelihood for each white dwarf is calculated 
using Equation 21 we use B ayes' Theorem as well as a 
Gaussian prior on the mass with a mean of 0.53 Mq and 
a standard deviation of 0.02 Mq, and a uniform prior on 



the temperature, to calculate the likelihood of the model 
parameters (m - M) Q and E(B-V) for each individual ob- 
ject, given the data. This step is described in Equation[5] 

Lj(((m-M) ,E(B- V)) | data) oc(5) 

p CO pOC 

I I Lj(data | (T cff ,(m- M) Q , mass, E(B - V))) 
Jo Jo 

(0.53 -mass) 2 ' 



*exp 



2(0. 02) 5 



dmass dT ff . 



We then combine the distance modulus-reddening like- 
lihood of all 59 WDs described in Equation [6] 



L(((m - M) , E(B - V)) | data) 



(6) 



J| J Lj(((m-M) ,E(B-V)) | data) 



and marginalise over the reddening parameter with a 
Gaussian prior with a mean of 0.04 and a a — 0.02. 
From this, we obtain the distribution function for the 
cluster as a whole, described in Equation [7] and shown in 
Figure El 



o 
o 




13.30 13.35 13.40 
(m - M) 



13.45 



Fig. 7. — The combined distribution function of distance modulus 
for the sample of 59 WDs. The peak of the distribution at 13.36 is 
the unreddened distance modulus value we quote for 47 Tuc along 
with its statistical uncertainty of 0.02. 



Our best fit unreddened distance modulus using these 
priors on mass and reddening is 13.36 ± 0.02. 



*exp 



L(((m - M) ) | data) = 
L((m - M) , E(B - V)) | data) 

(i 

(0.04-E(B-V)) 2_l 



(7) 



2(0.02) 2 



dE(B - V) 



Additionally, we consider how our distance estimate 
depends on our assumptions about the mass and redden- 
ing. To do this, we can consider our original four param- 
eter likelihood grid for each of our 59 white dwarfs. After 
marginalising out the effective temperature for each WD 
individually and combining the likelihood of all 59 WDs 
in Eqn. [51 we are left with the likelihood values on a grid 
over reddening, mass, and distance modulus. For each 
combination of reddening and mass on the grid, there is 



SEDs of White Dwarfs in 47 Tuc 



7 



a corresponding one-dimensional distance modulus likeli- 
hood distribution. The most likely distance modulus can 
then be found for each mass and reddening pair. The re- 
sult is a surface of most likely distance modulus values 
that is well approximately by a plane, described as 

(m - M) = 14.634 - (2.573 * mass) (8) 
+ (2.552 *E(B -V)) 

where the mass is in solar masses. This approximation 
to our likelihood grid is accurate to 0.003 mag and is 
valid only within the ranges of mass (0.50 — 0.56 M Q ) 
and reddening (0 — 0.08). 

5.2. Fitting the SEDs with DB Atmosphere Models 

Here we test the assumption that the WDs in our 
sample are DA atmosphere WDs by fitting our sam- 
ple of 59 WDs instead with DB atmosphere models 
i|Bergeron etalll2011l ). The SEDs fit with the DB models 
return temperatures that are on average 2200 K cooler 
than the fits with the DA models. These cooler temper- 
atures result in objects that are fainter and thus our true 
distance modulus is closer at a value of 12.90. If we com- 
pare the likelihoods of these fits, 53 of our 59 WDs have 
a higher likelihood value for our best fit DA model, with 
48 being better by a factor of 10 or more and 22 being 
better by a factor of 10 4 or more. The 6 remaining WDs 
that were fit better by a DB model have an average like- 
lihood that is better by a factor of 2. Additionally, the 
standard deviation of the differences between the best fit 
model and all 59 WDs is 0.11 mag for the DA models 
and 0.18 mag for the DB models. 

In t he field population of WDs, iTremblav fe Bergeronl 
(2008) found the ratio of DB to DA to be approximately 
25% in the temperature range of our WD sample. We 
test the effect on the determined distance modulus if 25% 
of our sample were instead DB WDs, by randomly fit- 
ting 25% of our WDs by DB models and the rest with DA 
models. Under these conditions, the distance modulus is 
13.32 with a random uncertainty of 0.06, easily encom- 
passing our distance modulus assuming all DA WDs. 

5.3. Error Analysis 

The statistical error of 0.02 that is quoted from our 
maximum likelihood analysis in Section 15.11 are lcr up- 
per and lower errors for the distance modulus calcu- 
lated directly from the summed likelihood distribution 
in Fig. [7J The systematic errors were determined en- 
tirely from the u ncertainties in our p hotometric calibra- 
tion taken from iKalirai et ahl ([201 ll ). We generate 10 
synthe tic WD SEDs from the models of [Trembla v et ahl 
(|20T1 evenly spaced between 8000 - 32000 K, spanning 
the range of effective temperature of our WD sample. We 
performed the same likelihood analysis as in Section [5. II 
using the uncertainties in the photometric calibration as 
the errors on each magnitude. For each of the 10 WDs, 
we extract a distance modulus from their likelihood dis- 
tributions along with the width of the distribution as the 
error. The widths of the likelihood distributions must be 
due solely to our photometric calibration uncertainties, 
as these synthetic objects come directly from the model 
grids. Figure [5] shows the systematic errors as a function 
of effective temperature for the synthetic WDs fit with a 
cubic spline function. For each of our 59 WDs, we used 
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Fig. 8. — The systematic uncertainties of the distance modulus 
as a function of effective temperature (T G ff ) for 10 synthetic WDs 
(triangles) fit with a cubic spline function. Inset: The distribution 
of systematic errors for the 59 WDs in our sample determined from 
the spline fit. Our quoted systematic uncertainty for the distance 
modulus of 47 Tuc is the mean of the sample of 59 WDs. 

the effective temperature to obtain the systematic error 
on its determined distance modulus from the spline fit 
in Fig. [8]For the final calibration error on the distance 
modulus of 47 Tuc, we quote the mean systematic error 
of the 59 WDs. 

Using this technique, we have determined the 
true distance modulus to 47 Tuc to be 13.36 ± 
0.02(random)±0.06(systematic) corresponding to a dis- 
tance of 4.70 ± 0.04(random)±0.13(systematic) kpc. An 
unknown fraction of DB objects in our sample will sys- 
tematically change this result, giving us a smaller dis- 
tance modulus, and wider error bars as the fraction in- 
creases. The closer distance results from the fact that 
DB white dwarfs of roughly the same SED shape are 
brighter than DAs, meaning we would have to fit them 
closer in order to match the observed magnitudes. The 
wider error bars result from the fact that DB models al- 
most universally fit our objects worse than DA models. 
Additionally, we don't know which objects are DBs, and 
so this causes a systematic error that needs to be ad- 
dressed by iterating our fitting routine while randomly 
assigning some fraction of our objects to be fit as DBs 
rather than DAs each time. The fit distance modulus and 
errors both change roughly linearly with DB fraction, as 
shown in Figure |H At 0% DBs we have the quoted fit 
values, at 25% DBs we have 13.32 ± 0.06 ± 0.06. 

We have tried a number of variations in order to as- 
sess how the distance modulus changes by varying some 
of our input assumptions. The white dwarfs forming 
in globular clu sters today should have masses between 
0.51- 0.55 (jRenzini fc Fusi Peccil9l^lRenzmi et al.l 
I1996T ). As discussed in Section l5T| we have used a Gaus- 
sian prior on the mass with a mean of 0.53 M Q and a a = 
0.02 in our distance determination. Based on Eqn[8j if we 
instead assume a WD mass of 0.50 M Q and 0.56 M Q , our 
unreddened distance moduli becomes 13.45 and 13.30, re- 
spectively. The reddening towards 47 Tuc has estimates 
ranging between E(B - V) = 0.024 dGratton et al.ll2003H 
and E(B - V) = 0.055 (jGratton et al.lll997f ) in the litera- 
ture. If we fix the reddening values to E(B - V) = 0.024 
and E(B - V) = 0.055, we obtain unreddened distance 
moduli of 13.33 and 13.41, respectively. Despite various 
changes to our initial assumptions, these are all consis- 
tent with our distance modulus of 13.36 given the random 
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Fig. 9. — The best fit distance modulus shown with varying frac- 
tions of DB/DA WDs in our sample. The uncertainties shown are 
the random and systemic uncertainties added in quadrature. 
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Fig. 10. — The distance modulus determined from our maximum 
likelihood technique shown as a function of reddening for various 
slices in WD mass space. 
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Fig. 11. — The previously determined unreddened distance mod- 
uli of 47 Tuc including the value from this study, plotted with their 
statistical uncertainties. The distances plotted from right to left are 
listed in Table IT1 from top to bottom. Our new measurement is the 
final point on the right. The different symbols represent the various 
methods used in the distance determination: main sequence fitting 
(open squares), horizontal branch fitting (filled circles), RR Lyrae 
(open circles), tip of the red giant branch (closed squares), methods 
using the white dwarfs (filled triangle), cluster kinematics (cross), 
and eclipsing binaries (open triangles). The solid line is the mean 
of the 22 distance moduli measurements (< (m - M) >= 13.30) 
and the dashed lines are one standard deviation (a = 0.13). 



uncertainties of 0.02 on each distance and the systematic 
uncertainty of 0.06 quoted above. We show in Figure ITOl 
how the distance modulus changes as a function of red- 
dening taken at different WD masses, all spanning the 
ranges in the literature. Again, we find all reasonable 
assumptions of the initial parameters to be consistent 
within our quoted uncertainties. 



As a summary, we show the previous distance deter- 
minations to 47 Tuc in Table Q] including the distance 
determined here in Figure 111! These values have been 
dereddened, where appropriate, using E(B - V) = 0.04 
and A v j E(B - V) = 3.1. We find the mean of these 
measurements to be < (m - M) >= 13.30 with a stan- 
dard error of the mean of 0.03 and a standard deviation 
of 0.13. 

6. DISCUSSION 

6.1. A Search for Infrared Excess 

Our dataset also provides us the opportunity to search 
for evidence of extrasolar planetary systems. A previ- 
ous searc h for planets in the c ore of 47 Tuc was per- 
formed bv lGilliland et al.l (|2000l ) using a time-series anal- 
ysis over 8.3 days searching for close in Jupiter-sized 
planets around ~ 34000 main sequence stars. They 
concluded that no light curve in their sample strongly 
suggested the presence of a planet, indicating that the 
planet frequency in 47 Tuc may be an order of magni- 
tude less than in our solar neighborhood. This may be 
expected in the dense stellar environment of the core of 
47 Tuc which could lead to orbit disruption and/or hinder 
planet survival. In addition, the planet-metallicity cor- 
relation indicates that the frequency o f planets increases 
with the host star's metallicity (e.g. iFischer fe ValentH 
|2005[) . making it even less likely to find a planet around 
a star with subsolar metallicities in 47 Tuc. 
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Fig. 12. — The observed minus the expected model magnitudes 
as a function of effective temperature (T e g ) for our WD sample in 
the four filters F390W (upper left), F606W (upper right), F110W 
(lower left), and F160W (lower right). 



Our dataset enables us to instead search for evi- 
dence of extrasolar planetary systems with the SEDs 
of our selected WDs. WDs have essentially metal- 
free a tmospheres (jPaauette e t al. 1986; Zuckcrma n et al.l 
120031 ) and low luminosities, allowing for the detection 
of cool low mass counterparts, such as debris or cir- 
cumstellar disks, as well as planets or brown dwarfs 
(lProbstlll983t iZuckerman fc Beck"hnlll987at TFarihi et al.l 
I2005t iBurleigh et al.l 1200^ 7" One form of evidence for 
the e xistence of these features is the presence of IR ex- 
cess (|Zuckerman fc Becklmlll987bt iGraham et all Il990t 
[ Kilic et all 120061 : Ivon Hippel et al.l 120071 : iFarihi et al.l 
120081 among others). The likelihood of a WD having 
an IR excess is correlat ed with its T e g a nd/or presence 
of photospheric metals (Fari hi et al.ll2009l ). the most well 
studied of which is calcium. The pollution of hydrogen- 



SEDs of White Dwarfs in 47 Tuc 



9 



rich photospheres of WDs with metals could be caused 
by circumstellar accretion. If a planetesimal or aster- 
oid is tidally dis r upted , a circumstellar dust disk could 
arise () Jural 120031 120081 ) and likely pollute the WD's at- 
mosphere with detectable metals. 

With no prior information about the calcium abun- 
dance of the WDs in our study, we examine their T e ff. 
Of our 59 WDs, all but 3 have fitted T eff > 10000 K, with 
the remaining three having T c ff = 9684 ± 1041 K, T off = 
9804 ± 921 K, and 9924 ± 841 K. We can estimate how 
many cases of IR excess we expect to see in our sample 
using upper and lower bounds provide d by various stud- 
ies in the literature. Recent work by IZuckcrma n et all 
(2010) show that ~ 30% of the photospheres of white 
dwarfs with T e ff between 13500 — 19500 K show evidence 
for being polluted with metals. In the same tempera- 
ture range, approximately 20% of DA WDs with metal 
bearing photospheres also show evide nce for warm cir- 
cumstellar material ()Kilic et al.l l2008f) . This indicates 
a fraction of dusty WDs to be 6%, which we take as 
an upper limit. A lower limit to th e fraction of warm 
DA WDs with dusty disks is - 1% dFarihi et al.l (2 009). 
obtained from recent lar ge WD surveys ([Liebert et al. 
2005tlKoester et al.1120051: iFarihi et al.ll200a iHoard et al.l 
20071 IMcCook fc Sionl 11999ft . From our 171 (59) WDs, 
we therefore expect between ^2 — 10 (^0 — 4) WDs to 
have IR excess. However, our reddest photometric mea- 
surement is F160W, a wide-H band filter with a central 
wavelength of 1536.9 nm. This is just where the spectral 
energy distribution of a typical warm WD would begin 
to show evidence for an IR excess given a companion or 
disk with T e fj ~ 800 K. In addition, we must also con- 
sider our photometric uncertainties. It thus seems un- 
likely that we would find clear evidence for an IR excess 
in our best sample of WDs. However, if we do find indi- 
cation for IR-excess in the H band, it could be confirmed 
with additional observations in redder filters. 

We have searched the SEDs in Figure |5] of all our WD 
candidates for evidence of an IR excess, focusing on the 
examination of the F160W band. Any evidence that we 
found for an enhanced F160W measurement appeared 
to be due to contamination by neighboring stars, diffrac- 
tion spikes, or cosmic rays. In addition, for each filter, we 
plot the observed magnitude minus the expected model 
magnitude as a function of effective temperature for each 
of our 59 WDs, shown in Figure [121 The few outlying 
points in the F160W filter at lower temperature are ob- 
jects that are too faint compared to the models, driven by 
large photometric uncertainties. If IR excess is present, 
we would expect to see points with negative values on the 
y-axis. Figure [T2l indicates that in the F160W band, our 
WDs are at most in excess of 0.2 magnitudes compared 
to the model. We question whether this could indicate 
a cool brown dwarf companion or debris disk. Using the 
evolutionary model s for cool brown dw arfs and exoso- 
lar giant planets of Baraffc et al. (2003), we combine a 
typical WD in our sample with the contribution from a 
10 Gyr dust free cool brown dwarf. Here, we convert the 
model IR data from the C IT system to that of 2MASS us- 
ing the transformations of lCarpen tcr (2001), as 2MASS is 
a nearly identical photometric system to that in WFC3 
and provides the closest comparison possible. In Fig- 
ure 1131 we show that for masses of a companion below 
0.06 M (~ 60 Mj), we would be unable to detect it. 
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Fig. 13.— The observed F606W, F110W, and F160W magni- 
tudes for a typical WD in our sample (solid squares). Also shown 
are the added contribution to these magnitudes from a 10 Gyr dust 
free cool brown dwarf companion with masses of 0.05 Mq (~ 50 
Mj) as open blue squares, 0.06 Mq (~ 60 Mj) as solid red tri- 
angles, 0.07 Mq (~ 70 Mj) as open green triangles, and 0.075 
Mq (~ 75 Mj) as solid cyan circles from the Baraffc et al. (2003) 
models. Here, we have converted the model mag n itudes to that of 
2MASS using the transformations in Carpenter (200l|). In order 
for us to detect a cool brown dwarf companion from our sample, 
the companion must be at least 0.06 Mq, which is an excess of 0.23 
magnitudes in the F160W filter. A color version of this figure is 
available in the electronic version. 

However, a cool brown dwarf with 0.06 Mq, may be de- 
tectable, as this leads to a difference in magnitude in 
the F160W band of > 0.23 magnitudes. Examining Fig- 
ure H2J it is clear that none of our 59 WDs have F160W 
in excess of 0.2 from the models. We conclude that we 
do not find any case in our sample that shows strong ev- 
idence for a genuine IR excess. Our search for planets 
in 47 Tuc continues, however, with this cycle 17 dataset. 
We have time series data of two fields outside the core 
of 47 Tuc, one is the stare field described in Section [5] 
consisting of F110W and F160W observations over ~ 20 
days with 10000 stars, and the other field has ACS images 
in F606W and F814W filters spanning ~ 9 months with 
20000 stars. Our analysis, currently underway, will be 
able to put a constraint on the binary frequency within 
47 Tuc and sear ch for exoplanets by a nalysis of the stel- 
lar light curves (|Woodlev et al.ll2012D . Not only will we 
examine the light curve of main sequence stars down to 
^0.2 Mq, but we will also be able to search for planets 
around the brighter WD population. A small radius ratio 
between the WDs and giant gas planets will significantly 
aid in our planetary search. 

6.2. Summary 

With this new HST dataset of the Galactic glob- 
ular cluster 47 Tuc, we have examined the SEDs 
of a sample of 59 WDs that are devoid of pho- 
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tometric contamination by neighboring stars or im- 
age defects. From this sample of SEDs, we have 
applied a novel approach to simultaneously estimate 
the distance of 47 Tuc from the comparison of 
theor etical model WD atmospheres (jTremblav et al.1 
1 2 1 If ) . We obtain a true distance modulus of 13.36 ± 
0.02(random)±0.06(systematic) corresponding to a dis- 
tance of 4.70 ± 0.04(random)±0.13(systematic) kpc. 
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TABLE 2 

Data on the Selected White Dwarfs 



ID 


Field 


R.A. 


Deal. 


F390W 


CF390W 


F606W 


CP606W 


F110W o-Fiiow 


F160W 


CP160W 


T cff a 


CT T eff 

(K) 






(J2000) 


(J2000) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) (mag) 


(mag) 


(mag) 


(K) 



WD 1 


Swath 1 


00 


:22 


:40 


.912 


-72 


:10: 


32 


.52 


25. 


.396 





.0351 


25. 


.362 





.0320 


25 


.500 





.1301 


25. 


.515 


0. 


.2238 


12047 


1602 


WD 2 


Swath 1 


00 


:22 


:25 


.686 


-72 


:10 


11 


.84 


23. 


.620 





.0207 


23 


.950 





.0194 


24 


.473 





0464 


24. 


.430 


0. 


0705 


20698 


1361 


WD 3 


Swath 2 


00 


:21: 


:57 


.095 


-72 


:10: 


18. 


59 


24 


263 





.0223 


24 


,381 





.0213 


24 


.866 





.0623 


24. 


.898 


0. 


1030 


16132 


1041 


WD 4 


Swath 2 


00 


:22 


:00 


.123 


-72 


:10: 


12 


.02 


24. 


.399 





.0163 


24 


,495 





.0182 


24 


.722 


0. 


0566 


24 


.732 


0. 


0987 


14570 


720 


WD 5 


Swath 2 


00: 


:22 


:00 


.719 


-72 


:09: 


15. 


.09 


25. 


.217 





.0270 


25 


125 





.0322 


25 


.303 





.0763 


25 


.297 





1365 


11687 


1161 


WD 6 


Swath 2 


00 


:22 


:26 


.133 


-72 


:09: 


20 


01 


24 


.108 





.0134 


24 


.287 





.0235 


24 


.633 





0612 


24. 


631 


0. 


1141 


16372 


921 


WD 7 


Swath 2 


00 


:22 


:01 


.506 


-72 


:08: 


34. 


31 


23 


.314 





.0220 


23 


.684 





,0153 


24 


.070 





0386 


24. 


.062 


0. 


0683 


20458 


1281 


WD 8 


Swath 2 


00 


:22 


:25 


.423 


-72 


:09: 


09 


.83 


25. 


.257 





.0338 


25 


.138 





.0316 


25 


.325 


0. 


0828 


25. 


373 


0. 


1618 


11566 


1281 


WD 9 


Swath 2 


00: 


:22 


:23 


.633 


-72 


:09: 


00. 


06 


23. 


.728 





.0242 


24 


.001 





.0305 


24 


.552 





.0550 


24. 


.795 


0. 


1454 


20257 


1722 


WD 10 


Swath 3 


00 


:21: 


:31 


.487 


-72 


:08: 


35. 


91 


23. 


.950 





.0213 


24 


129 





.0239 


24 


.624 





0618 


24. 


.996 


0. 


1379 


17574 


1201 


WD 11 


Swath 3 


00 


:21: 


:39 


.401 


-72 


:09: 


03. 


.65 


25 


.658 





.0271 


25 


.498 





.0292 


25 


.538 





1380 


25. 


.140 


0. 


2300 


9684 


1041 


WD 12 


Swath 3 


00 


:21 


:51 


.751 


-72 


:09 


31 


15 


24, 


.491 





.0146 


24 


.527 





.0143 


24 


.726 





0558 


24. 


953 


0. 


1187 


13889 


640 


WD 13 


Swath 3 


00 


:22 


:01 


.247 


-72 


:09: 


16. 


43 


25. 


.190 





.0318 


25 


.088 





.0218 


25 


.168 





.0731 


25 


.056 


0. 


1274 


10966 


961 


WD 14 


Swath 3 


00 


:21 


:47 


.243 


-72 


:08. 


22. 


.44 


25 


.310 





.0232 


25 


.168 





.0257 


25 


.292 





.0747 


24. 


.953 


0. 


0945 


10285 


640 


WD 15 


Swath 3 


00 


:21 


:47 


.989 


-72 


:08. 


15. 


.89 


25 


613 





.0291 


25 


.463 





.0316 


25 


.542 





0982 


25. 


.239 


0. 


1113 


10045 


801 


WD 16 


Swath 3 


00 


:21: 


:41 


.681 


-72 


:07: 


:49. 


43 


24. 


.062 





.0179 


24 


.265 





.0182 


24 


.778 





0581 


24 


.658 


0. 


0921 


17374 


1001 


WD 17 


Swath 3 


00 


:21 


:42 


.075 


-72 


:07: 


47. 


.08 


24. 


.640 





.0263 


24 


.711 





.0233 


25 


.057 





0772 


25. 


144 


0. 


1618 


14730 


1041 


WD 18 


Swath 3 


00 


:21: 


:46 


.162 


-72 


:07: 


:54. 


27 


23. 


.548 





.0210 


23 


.877 





.0178 


24 


.352 





.0638 


24. 


.332 


0. 


.0914 


19977 


1401 


WD 19 


Swath 3 


00 


:22 


:02 


.034 


-72 


:08. 


'35 


65 


23. 


.332 





.0215 


23 


.688 





.0191 


24 


.209 





0448 


24. 


.371 


0. 


0739 


22220 


1481 


WD 20 


Swath 3 


00 


:21: 


:53 


.267 


-72 


:08. 


02. 


.38 


24 


.810 





.0203 


24 


.744 





.0225 


25 


.257 





0831 


25. 


.381 


0. 


2374 


13048 


961 


WD 21 


Swath 4 


00: 


:21: 


:30 


.986 


-72 


:08: 


00. 


.98 


24 


.132 





.0257 


24. 


.272 





.0247 


24 


.622 





.0529 


24. 


.548 


0. 


0971 


15812 


1041 


WD 22 


Swath 5 


00 


:21: 


:07 


.243 


-72 


:04: 


31. 


.54 


24. 


833 





.0275 


24 


.712 





.0213 


25 


.087 





0823 


24. 


.955 


0. 


1090 


12207 


1001 


WD 23 


Swath 5 


00 


:21 


:15 


.803 


-72 


:06 


00 


.35 


24, 


.392 





.0309 


24 


.438 





.0261 


24 


.956 





0803 


24. 


806 


0. 


1638 


14971 


1201 


WD 24 


Swath 5 


00: 


:21: 


XI. 


.805 


-72 


:04: 


:23. 


.82 


23. 


.039 





.0199 


23 


.469 





.0151 


24 


,030 





.0558 


23. 


.951 


0. 


0668 


23221 


1401 


WD 25 


Swath 5 


00 


:21 


:23 


.488 


-72 


:03: 


:44. 


.75 


23 


.803 





.0237 


24 


.057 





.0245 


24 


.460 





0606 


24. 


.654 


0. 


0875 


18816 


1361 


WD 26 


Swath 5 


00 


:21: 


:32 


.257 


-72 


:05: 


38. 


40 


22 


316 





.0209 


22 


.940 





.0134 


23 


.570 





0345 


23. 


.641 


0. 


0564 


31712 


2563 


WD 27 


Swath 6 


00: 


:21: 


:26 


.729 


-72 


:01: 


40. 


22 


25. 


.413 





.0376 


25 


.337 





.0198 


25 


.266 


0. 


.0843 


25 


.044 





1158 


9924 


841 


WD 28 


Swath 6 


00 


:21 


:23 


.115 


-72 


:03: 


■44. 


39 


23 


.816 





.0156 


24. 


.088 





.0360 


24 


.429 





0449 


24. 


.652 


0. 


0810 


18655 


1121 


WD 29 


Swath 7 


00 


:21 


:39 


.901 


-71 


:59: 


46 


.59 


24 


.809 





.0260 


24 


.789 





.0226 


25 


.096 





0698 


25. 


.418 


0. 


1630 


13569 


1041 


WD 30 


Swath 7 


00: 


:21: 


:31 


.482 


-72 


:01: 


:05. 


.77 


24, 


.259 





.0117 


24 


.505 





.0125 


24 


.822 





.0581 


25 


.031 


0. 


0938 


17614 


720 


WD 31 


Swath 7 


00 


:21 


:48 


.558 


-72 


:00 


42 


.08 


23 


.079 





.0356 


23 


.569 





.0165 


24 


.071 





0388 


24. 


.183 


0. 


0619 


25544 


2242 


WD 32 


Swath 8 


00 


:21 


:55 


.432 


-71 


:58: 


34 


.17 


24. 


.546 





.0294 


24. 


.549 





.0256 


24 


.842 





0576 


24. 


.642 


0. 


0772 


13329 


1041 


WD 33 


Swath 8 


00 


:22 


:00 


.315 


-71 


:58: 


30. 


.80 


25. 


.161 





.0322 


25. 


.080 





.0281 


25 


.279 





.0724 


24. 


.927 


0. 


0880 


11086 


961 


WD 34 


Swath 8 


00 


:21 


:59 


.219 


-71 


:58: 


.45. 


.40 


25 


.135 





.0303 


25 


.140 





.0240 


25 


.094 





0649 


25 


.171 


0. 


1203 


11246 


1121 


WD 35 


Swath 8 


00 


:21: 


:57 


.370 


-71 


:59: 


11 


32 


22 


.845 





.0233 


23 


.291 





.0186 


23 


.843 





0332 


23 


.997 


0. 


0509 


25584 


1642 


WD 36 


Swath 8 


00 


:22. 


:05 


.843 


-71 


:59: 


28. 


07 


24, 


.740 





.0295 


24. 


.749 





.0277 


24 


.944 





.0686 


24. 


924 


0. 


1045 


13209 


1161 


WD 37 


Swath 9 


00 


:22 


:21 


111 


-71 


:57. 


49 


36 


23 


.525 





.0305 


23 


.832 





.0247 


24 


.309 





0454 


24. 


645 


0. 


0911 


21419 


1802 


WD 38 


Swath 9 


00 


:22 


12 


.785 


-71 


:58: 


:54. 


71 


25. 


.127 





.0260 


25 


,015 





.0220 


25 


.241 





.0781 


25 


.100 


0. 


1078 


11446 


961 


WD 39 


Swath 10 


00: 


:22 


:55 


.003 


-71 


:57: 


:27. 


72 


24. 


.810 





.0188 


24 


.749 





.0178 


25 


.060 





.0700 


25. 


.658 


0. 


1799 


12888 


881 


WD 40 


Swath 10 


00 


:22 


:59 


.788 


-71 


:57. 


41. 


.86 


25 


.047 





.0274 


24. 


.976 





.0250 


25 


.020 





.0647 


24. 


.961 


0. 


1094 


11006 


881 


WD 41 


Swath 10 


00 


:22 


:36 


.414 


-71 


:57: 


.17. 


.74 


24. 


.791 





.0306 


24 


.775 





.0284 


24 


.988 





0719 


25 


.118 


0. 


.1012 


13409 


1161 


WD 42 


Swath 10 


00 


:22 


:56 


.251 


-71 


:57: 


55 


.80 


22. 


.920 





.0167 


23 


.420 





.0225 


23 


.819 


0. 


.0340 


23. 


777 


0. 


0557 


23461 


1201 


WD 43 


Swath 10 


00: 


:22 


:33 


.763 


-71 


:58: 


55. 


.87 


25. 


.673 





.0391 


25 


,548 





.0274 


25 


,462 





.0908 


25. 


.290 


0. 


1584 


9804 


921 


WD 44 


Swath 11 


00 


:22 


:59 


.738 


-71 


:57: 


,12. 


01 


25. 


.090 





.0386 


25 


,060 





.0247 


25 


.160 





0726 


25 


.228 


0. 


.2013 


11927 


1441 


WD 45 


Swath 12 


00 


:23 


:33 


.344 


-71 


:58. 


26 


.77 


24, 


.259 





.0177 


24 


.387 





.0193 


24 


.446 





.0725 


24. 


.447 


0. 


1015 


14690 


801 


WD 46 


Swath 12 


00 


:23 


:26 


.967 


-71 


:58: 


'57. 


.50 


23 


.188 





.0194 


23 


.617 





.0193 


24 


.063 





0596 


24. 


366 


0. 


1228 


23141 


1561 


WD 47 


Swath 12 


00 


:23: 


:32 


.372 


-71 


:59: 


:45. 


.70 


22 


.430 





.0133 


23 


.026 





.0158 


23 


.507 





.0376 


23. 


.823 


0. 


0538 


29509 


1762 


WD 48 


Swath 12 


00 


:23: 


:21 


.810 


-71 


:59: 


18. 


03 


24 


.233 





.0179 


24. 


.376 





.0250 


24 


.915 





.0882 


25. 


.064 


0. 


1247 


16653 


1121 


WD 49 


Swath 12 


00 


:23: 


:37. 


.390 


-72 


:00 


36 


89 


23. 


.387 





.0151 


23 


.753 





.0208 


24 


.122 





1033 


24. 


.002 


0. 


1056 


20017 


1361 


WD 50 


Swath 12 


00 


:23: 


:28 


.612 


-72 


:00: 


13 


.42 


24. 


.186 





.0227 


24 


.289 





.0250 


24 


.459 





0823 


24. 


.690 


0. 


1327 


14770 


1001 


WD 51 


Stare 


00: 


23 


:32. 


.099 


-72 


:09: 


39. 


30 


23. 


.982 


0. 


0346 


24 


.166 





.0248 


24 


.569 





.0471 


24. 


.632 


0. 


2060 


17254 


1481 


WD 52 


Stare 


00: 


:23 


:00 


.862 


-72 


:10: 


38. 


40 


25. 


293 


0. 


.0497 


25 


.221 





.0258 


25 


.361 





.0681 


25. 


.014 


0. 


0842 


10725 


921 


WD 53 


Stare 


00: 


23: 


:17. 


641 


-72 


:09: 


53. 


81 


25 


.215 


0. 


0370 


25 


.103 





.0244 


25 


.282 





0706 


25. 


.734 


0. 


2999 


11687 


1281 


WD 54 


Stare 


00: 


23 


:09 


.435 


-72 


:10: 


11 


30 


23. 


998 


0. 


0178 


24 


.209 





.0172 


24 


.527 


0. 


.0241 


24. 


668 


0. 


0740 


16973 


760 


WD 55 


Stare 


00 


:23 


:18 


.825 


-72 


:09: 


08. 


07 


23. 


.725 


0. 


.0201 


24. 


.025 





.0175 


24 


.597 





.0727 


24. 


.505 


0. 


0780 


19416 


1241 


WD 56 


Stare 


00 


:23 


12 


136 


-72 


:09: 


23 


12 


25 


264 





.0520 


25. 


.184 





.0252 


25 


.310 





0333 


25. 


.692 


0. 


1589 


12047 


1081 


WD 57 


Stare 


00 


23 


:03 


033 


-72 


:09: 


05 


57 


25 


514 


0. 


0426 


25 


467 





.0291 


25 


.430 





.0524 


25. 


.235 


0. 


1018 


10325 


720 


WD 58 


Stare 


00: 


:23 


13 


.822 


-72 


:08: 


18. 


48 


23. 


385 


0. 


.0147 


23 


.794 





.0198 


24 


.341 





.0291 


24. 


.294 


0. 


0418 


23341 


1001 


WD 59 


Stare 


00: 


:22 


:55. 


.764 


-72 


:08: 


56 


88 


22 


639 


0. 


0105 


23 


.219 





.0159 


23 


.758 





0129 


23. 


890 


0. 


0262 


29028 


801 



a Effective temperature estimates assume a white dwarf mass of 0.53 M©. 
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Fig. Set 3. Images of our WD sample 

Fig. Set 6. Spectral Energy Distributions and likelihood fits for our WDs 



